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Abstract 

We study the evolution of homogeneous and isotropic, fiat cosmological models within 
the general scalar-tensor theory of gravity with arbitrary coupling function and potential 
and scrutinize its limit to general relativity. Using the methods of dynamical systems for 
the decoupled equation of the Jordan frame scalar field we find the fixed points of flows 
in two cases: potential domination and matter domination. We present the conditions 
on the mathematical form of the coupling function and potential which determine the 
nature of the fixed points (attractor or other). There are two types of fixed points, both 
are characterized by cosmological evolution mimicking general relativity, but only one of 
the types is compatible with the Solar System PPN constraints. 
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1 Introduction 



Scalar-tensor theories (STT) of gravitation [TJ [2, [3] have emerged in different contexts of 
theoretical physics, e.g. in Kaluza-Klein type unifications, supergravity, and low energy ap- 
proximations of string theories. In cosmology STT has been invoked to model the accelerated 
expansion of inflation and dark energy. However, observations in the Solar System tend to 
indicate that in the intermediate-range distances the present Universe around us is successfully 
described by the Einstein tensorial gravity alone [H El E] . This means that only such models 
of scalar-tensor gravity are viable which in their late time cosmological evolution imply local 
observational consequences very close to those of Einstein's general relativity (GR) [TJ. 

The methods of dynamical systems provide natural tools to analyze the problem. In this 
paper, summarizing our recent work |8J, we take the Jordan frame and consider general scalar- 
tensor theories which contain two functional degrees of freedom, the coupling function u(^f) and 
the scalar potential V(^f). We perform the dynamical systems analysis for the flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) backgrounds with ideal barotropic fluid matter. Our 
strategy is to find the fixed points for the scalar field dynamics and compare these with 
the conditions of the limit of general relativity in the Solar System, as established by the 
parameterized post-Newtonian (PPN) formalism. Therefore, if the functional forms of u(ty) 
and V(*f>) are specified from some considerations (e.g. the compactification manifold), our 
results allow to determine the fixed points along with their type and thus immediately decide 
whether general relativity is an attractor, i.e. whether the model at hand is viable or not. 

2 Scalar-tensor cosmology as a dynamical system and 
the limit of general relativity 

We consider a general scalar-tensor theory in the Jordan frame given by the action functional 



S= 2^' d X ^ 
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Here ufi?) is a coupling function and V(\&) is a scalar potential, V M denotes the covariant 
derivative with respect to the metric g pv and S m is the matter part of the action as all other 
fields are included in Xm- In order to keep the effective Newtonian gravitational constant 
positive [9] we assume that < \P < oo. 

The field equations for the flat (k = 0) FLRW line element and perfect barotropic fluid 
matter, p = wp, read 
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Here H = a/a and we assume p > 0. Eqs. (J2J) — (EI) are too cumbersome to be solved analyti- 
cally, but useful information about the general characteristics of solutions can be obtained by 
rewriting (j2])-([5]) in the form of a dynamical system and finding the fixed points which describe 
the asymptotic behaviour of solutions. 

The phase space of the system is spanned by four variables H, p}. Defining ^ = 

x, ^ = y the dynamical system corresponding to equations (J2])-([2D can be written as follows: 
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Based on these equations we may make a couple of quick qualitative observations about 
some general features of the solutions. For example, the limit \P — ► in general implies 
\H\ — > oo, hence the solutions can not safely pass from positive to negative values of \I/ 
(from "attractive" to "repulsive" gravity), but hit a space-time singularity as the curvature 
invariants diverge. Similarly, the limit 2u + 3 — > implies \H\ — > oo with the same conclusion 
that passing through ufi?) = — | would entail a space-time singularity and is impossible. 
(Let us remark here, that these observations are quite general and do not preclude specially 
fine-tuned solutions in some fine-tuned models which may remain regular while crossing these 
points pig.) 



The limit 
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deserves a more closer examination. Let us define by (2a;(x*)+3) 1 



0. Expressing H from the Friedmann constraint 
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makes clear that \H\ diverg 
the latter case is determined by the first term under the square root above. We can compute 
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its limit by Taylor expanding 



lim (2u(x) + 3)y 2 



lim 7 

dx V2oj(x)+3 



" 2 sin 2 6 1 



1 



r cos 







tan 2 6> 



oo , 



if 
if 
if 



_d_ 

dx 



1 



2u;(x)+3 
1 

dx V2^(x)+3 



<!x Gu^xj+s) 



17 ' 2 dx 2 V2w(x)+3 
7^0, 

x=x* 

= 

x=x* 

= 



r 2 cos 2 # + 



<Lr2 (2a;(x)+s) 



dx 2 



(2w(x)+s) 



7^0 

= 



'11^ 



where Ax = rcosO, Ay = rsin6* was taken. (We have neglected the unphysical direction 



16*1 = | that corresponds to approaching the point 
where \H\ is divergent.) So, in the limit x 
lowest non-zero derivative of 1 



x*,y = 0) along the line x 



the value of H is determined by the 
If the both the first and second derivative vanish, then 
\H\ diverges implying a spacetime singularity. If the first derivative vanishes but the second 
derivative is not zero, 
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then H is finite but (possibly) different for each solution as it depends on the angle of approach 
6, while the Friedmann equation in this case acquires an extra term when compared to general 
relativity. If the first derivative is not zero, 
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then H approaches the value if 2 = £r(p + V(x ir )), mimicking the Friedmann equation of 

2 2 

general relativity with 8nG = — and A = —V(x+). 

To summarize, we have just observed that in the limit (a) 2 u(x)+3 ~ > ^> O 3 ) V ~ * ^ ^ ne 
Friedmann constraint ([TO]) tends to the form of general relativity if (c) ^ 2 w(x 1 )+3) 2 ^x ^ ^- 
It must be also emphasized here that the process of taking the Taylor expansion (fTTj) hinges 
on the assumption that (d) 2c ^ +3 is differentiable (derivatives do not diverge) at x*. In this 
context one may also ask when the full set of equations ([6])- (|9j) attains the form of general 
relativity. It is easy to see that besides (a)-(d) one must also impose 
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but the latter is automatically satisfied if (c) holds, due to (fill) . (|T3|) . Therefore we may 
tentatively call the conditions (a)-(d) 'the general relativity limit of scalar-tensor flat FLRW 
cosmology'. 
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It is interesting to compare the cosmological GR limit to the GR limit obtained from PPN, 
which characterizes the slow motion approximation in a centrally symmetric gravitational field. 
Although the mathematical assumptions underlying the PPN formalism are clearly different 
from our cosmological reasoning above, we may still ask whether the results of both schemes 
agree with each other. In the context of PPN it is well established that the solutions of 
scalar-tensor theory approach those of general relativity when [9] 

1 o, ,„ ,:,„„ ^~o. 



2w(x) + 3 ' (2u(x) + 3) 3 dx 
Comparison shows that the cosmological conditions (a)-(d) are marginally stricter than the 
PPN condition (|T5l) . since (a), (c), (d) imply that (|T5l) is satisfied, but ( TT5l) does not necessarily 
guarantee that (c) or (d) holds. 

Let us also note that there is also another special case x,, realized at 

dV(x) 

p = 0, y = 0, 2V(x.)-x. — ^ =0, (16) 

dx 

when the cosmological equations (jEJ)-© relax to those of general relativity featuring de Sitter 
evolution. However, as the value of w(x,) is not fixed by the condition ([TBI) , this case does 
not conform with the GR limit of PPN. Therefore, even when the limits (a)-(d) and ([TBI) can 
be cosmologically indistinguishable, Solar System observations in the PPN framework can in 
principle reveal which of the two is actually realized. (In this paper when using the phrase 
'the GR limit of STT we mean the conditions (a)-(d), as these take the STT cosmological 
equations to those of general relativity and also guarantee that the PPN condition is satisfied. 
But note that some authors pi], [12] have not necessarily used the same definition.) 

The general relativity limit of STT is purely given in terms of x and y. In the following 
we extract from the full dynamical system ([6])- (J9j) an independent subsystem for {x,y}, find 
its fixed points and check whether the limit of general relativity matches to an attractive fixed 
point. 

3 Fixed points for potential domination (V ^ 0, p = 0) 

In the four phase space dimensions = x, ^ = y,H,p} the physical trajectories (orbits 
of solutions) are those which satisfy the Friedmann constraint (JSJ). In the limit of vanishing 
matter density the phase space shrinks to three dimensions {x = ty,y = where the 

Friedmann constraint restrains the physical trajectories to span two dimensions. We may 
solve the Friedmann constraint for H, as ( [T0[) . substitute it into Eq. and thus in effect 
reduce the system 2-dimensional: 

x = y 



y = (i - Mkr-s t) y 2 { +) h^M^) + % 2 + i2«w(x) y (17) 
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Condition 



Eigenvalues 
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Table 1: Fixed points and their eigenvalues for the V ^ 0, p = case. 

This constitutes a projection of the trajectories on the original two-dimensional constraint 
surface in (x,y,H) to the (x,y) plane. The projection yields two "sheets": the "upper sheet" 
marked by the ~ sign, and the "lower sheet" marked by the the ^_|_^ sign in Eq. (|17|) . 

Standard procedure reveals that the dynamical system (ITTj) is endowed with two fixed 
points, Table [1] lists their conditions and eigenvalues. The first fixed point satisfies [131 E] 

^| # ».-m«.)=0, (18) 

which matches the second limit (fIBj) . discussed in Sec. 2. The second fixed point ^ satisfies 

1 1 dcu 







2cj(^) + 3 (M**) + 3) 2 d^ 



i.e. exactly the same conditions (a)-(d) as the limit of general relativity for flat FLRW STT 
cosmology, discussed in the end of Sec. 2. 

From Eq. (TTOT) it is straightforward to compute that the values of H corresponding to the 
fixed points and \P* are H, = r—\y~^y~ an d H+ = \J ^-^7^, respectively. The 
result, which mimics de Sitter evolution in general relativity, was expected, since we saw in 
Sec. 2 that under the fixed point conditions (ITS!) and (IT31) the full STT equations ©-(JED reduce 
to the equations of general relativity. 

4 Fixed points for matter domination (V = 0, p ^ 0) 

In the case of cosmological matter (p > 0) and vanishing scalar potential the Friedmann con- 
straint restricts the solutions onto a three-dimensional surface in four phase space dimensions 
(x = = ^, H, p). However, the system is amenable to a change of the time variable [7j that 
allows to combine the field equations into a dynamical equation for the scalar field which does 
not manifestly contain the scale factor or matter density. In the Jordan frame this amounts 
to defining a new time variable, dp = h c dt = H + ^ dt , and deriving from Eqs. ([2])-(j4|) 
a "master" equation for the scalar field [151 IS] • It can be written in a form of a dynamical 
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system for variables \I/ = x, ^f' = z (f = df /dp) 



x' 
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The signs in Eq. ( 120]) correspond to the "upper" and the "lower" sheet as before in accordance 
with the constraint equation (JTU1) . 

The "master" equation can be relied on as long as h c = |i? + 7^7 1 is finite. At \I/ = the 
quantity h c diverges making the t-time to stop with respect to the p-time. Hence all t-time 
trajectories with finite \& get mapped to = 0, giving a false impression of a fixed point 
there. However, in Sec. 2 we concluded that — comes with a space-time singularity and 
exclude it from present analysis. 

The other problematic points can be discussed by noticing that in terms of the new time 
variable p the Friedmann constraint ([2]) can be written as 

h 2 = ^ (21) 

C 3* (l - ^1±^ X ' 

To keep h c real, the right hand side of Eq. fT2~T]) must be nonnegative, thus constraining the 
dynamically allowed regions of the two-dimensional phase space of the scalar field. It 

is also important to verify that a fixed point in the p-time does indeed correspond to a fixed 
point in the cosmological t-time. But since on physical grounds it is reasonable to consider only 
the trajectories with finite h c , it is immediate that ty' = implies \l/ = due to redefinition 
p = p(t). 

Let us consider the cosmological matter behaving like dust (w = 0). An argument com- 
pletely analogous to the one put forth for Eq. (fTTl) . reveals a single fixed point, satisfying 
the conditions (a)-(d) dubbed as the limit of general relativity for flat FLRW STT. The cor- 
responding eigenvalues, to be evaluated at the fixed point coordinate, are given in Table [2j 
In particular, this point is an attractor on the "upper" sheet if 4^ > 0, while on the "lower" 
sheet attractor behavior is not possible. From Eq. ffTOl now it also follows that at the fixed 
point the evolution of the universe obeys the usual Friedmann equation from general relativ- 
ity, Hi = This is expected, as the fixed point conditions were identical to the general 
relativity limit. 

The dynamical system in the radiation dominated regime (w = |) has no fixed points. For 
small values of the system is ruled by friction on the "upper" sheet, as the ~ sign of the 
dominating term forces the vector flow to converge to the ty" = axis. On the "lower" sheet, 
the the effect is the opposite (anti- friction). 
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Case 



Fixed point 



Condition 



Eigenvalues 
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Table 2: Fixed points and their eigenvalues for the p 0, V = case. 



5 Conclusion 

We have considered flat FLRW cosmological models in general scalar-tensor theories with 
arbitrary coupling function oufi?) and scalar potential V(^ r ) in the Jordan frame. Using the 
methods of dynamical systems we have found the scalar field fixed points in two distinct 
asymptotic regimes: potential domination (V ^ 0,p = 0), and matter domination (V = 
0, p ^ 0). In nutshell there are two types of fixed points arising from different mechanisms: 

from a condition on the potential and from the singularity of the scalar field kinetic 
term. Approaching both types of fixed points the cosmological equations coincide with those of 
general relativity, yielding de Sitter expansion in the potential domination case and Friedmann 
evolution in the matter domination case. However, for the Solar System experiments in the 
PPN framework only the fixed points of type give predictions identical with those of general 
relativity. The nature of fixed points (attractor or otherwise) depends on the functional forms 
of u(^f) and V(^) according to corresponding eigenvalues given in Tables 1 and 2. Therefore, 
in Jordan frame analysis, general relativity is an attractor for a large class of scalar-tensor 
models, but not for all. 

Provided the transformation relating the Jordan and the Einstein frame is regular, there 
is an exact correspondence between the two frames and the Jordan frame phase space results 
should carry over to the Einstein frame [T7J [HI EE 120] • This is the case for the fixed point 
ty,. However, as the transformation of the scalar field fails to be regular in the limit of general 
relativity, the properties of the ^ fixed point may be altered in the Einstein frame [TTj . To 
establish whether or how the correspondence holds in this case calls for a separate matching 
investigation in the Einstein frame. 
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